Unified paradigm for interface dynamics 
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In this paper we develop a common theoretical framework for the dynamics of thin featureless 
interfaces. We explicitly demonstrate that the same phase field and velocity dependent one-scale 
models characterizing the dynamics of relativistic domain walls, in a cosmological context, can also 
successfully describe, in a friction dominated regime, the dynamics of nonrelativistic interfaces in 
a wide variety of material systems. We further show that a statistical version of von Neumann's 
law applies in the case of scaling relativistic interface networks, implying that, although relativistic 
and nonrelativistic interfaces have very different dynamics, a single simulation snapshot is not able 
to clearly distinguish the two regimes. We highlight that crucial information is contained in the 
probability distribution function for the number of edges of domains bounded by the interface 
network and explain why laboratory tests with nonrelativistic interfaces can be used to rule out 
cosmological domain walls as a significant dark energy source. 



I. INTRODUCTION 

Phase field models (PFM) have been applied with suc- 
cess to describe the interface dynamics of various nonrel- 
ativistic systems, including foams [TJ [5] and grain growth 
[3] (see also H E] for recent reviews of PFM). In PFM 
sharp interfaces are modeled as diffuse boundaries, with 
a thickness much smaller than the curvature radius, sep- 
arating regions characterized by different values of the 
order parameters. In the simplest cases with noncon- 
served order parameters, such as ideal grain growth, cel- 
lular patterns bounded by an interface network tend to 
evolve toward configurations with less interfacial energy 
(or equivalently, less surface area). This leads to the 
gradual elimination of entire domains, resulting in the 
increase of the characteristic scale L with physical time t 
as L (X t^^'^ both for two- and three-dimensional networks 

In these systems the velocity v is proportional to the 
mean curvature of the interface at each point. This is 
the crucial ingredient behind von Neumann's law 8J for 
the evolution of the area of individual domains in two- 
dimensional interface networks with (free) Y-type junc- 
tions: dttn/dt (X (n — 6), where a„ is the area of a domain 
with n edges. Since the average number of edges of an in- 
terface network with Y-type junctions is 6, the only pos- 
sible static configuration consistent with von Neumann's 
law is a hexagonal lattice. Note that if the interfaces 
have a finite width, von Neumann's law may not apply 
to individual domains but is reasonably obeyed on aver- 
age [S]. Extended versions of von Neumann's law exist 
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in both two [10] and three dimensions [TTI - fT5] (see also 
Ref. [H]). Although other laws, such as Lewis's law and 
Aboav's law, have also been used in the characterization 
of classical patterns [TIHTB] , in this paper we shall focus 
on von Neumann's law, which is more directly relevant 
for the evolution of an interface network. 

On the other hand, the dynamics of relativistic inter- 
face networks, considering various possible types of de- 
fect junctions, has been studied in detail in a cosmolog- 
ical context, where domain walls have been proposed in 
the past as an interesting dark energy candidate [T7]. In 
Refs. im US] , the authors have used large and accurate 
sets of field theory numerical simulations in combination 
with semi-analytical analysis to provide very strong ev- 
idence for a no-frustration conjecture, which invalidates 
domain walls as a viable dark energy candidate. 

The main aim of this paper is the development of a 
unified theoretical paradigm for interface dynamics which 
includes both relativistic and nonrelativistic systems in a 
unified framework. We set up a velocity-dependent one- 
scale (VOS) model and use it to describe the evolution of 
the characteristic length L and velocity v of thin, feature- 
less interface networks both in the relativistic and nonrel- 
ativistic limits. We show that the same models which de- 
scribe the dynamics of cosmic relativistic interfaces also 
account, in a friction dominated regime, for the dynam- 
ics of cellular patterns in many material systems, such as 
soap froths, lipid monolayers, and metals. Then we com- 
pare single simulation snapshots of high resolution rela- 
tivistic and nonrelativistic interface network simulations, 
discussing the cosmological implications of our findings, 
in particular for dark energy. We also demonstrate that 
von Neumann's law still applies in the relativistic regime, 
albeit only on average. In this paper we use natural units 
with c = h = 1 and the Einstein summation convention. 
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II. UNIFIED FRAMEWORK 

We shall consider a multicomponent field theory La- 
grangian with a discrete symmetry. For simplicity let us 
choose a model with a scalar field multiplet (j^e described 
by the Lagrangian 



C = X-V{^e): 



(1) 



where X = — 9^(/)e 9^0e/2 and V{(j)e) is the potential 
(with the exception of Eq. ( f 3 ), in the following we shall 
omit the index e). This model admits interface solu- 
tions if the potential V{(l)) has at least two discrete de- 
generate minima. Although models described by more 
general Lagrangians can be considered, that does not in- 
troduce additional dynamical degrees of freedom, as long 
as the interfaces are thin and featureless [ini [20] ■ Still, 
one should bear in mind that there are very rich spectra 
of geometrical, energetic, and topological network prop- 
erties associated with different choices of the potential 

Varying the action. 



S = 



(2) 



with respect to the scalar field (p one obtains the following 
equation of motion 



dV 



(3) 



Here g = (let{ga[j) and g^p is the metric tensor. In a 
flat homogeneous and isotropic three-dimensional back- 
ground the line element is ds^ = —dt^ + i?^(t)dx • dx, 
where R{t) is the scale factor, t is the physical time, and 
X ~ (x, y, z) are comoving coordinates. In this case, the 
equation of motion for the scalar field (f), given by Eq. 
([3]), can be written as 



dV 



(4) 



where a dot represents a derivative with respect to phys- 
ical time, l^^ ^3H, H ^ R/R, r = i?(t)x, and is the 
physical Laplacian. In the presence of frictional forces the 
characteristic damping length becomes = 3H + i^^: 
where £{ is the friction length scale. 

It has been shown in Ref. 20 that if the thickness 
of the interface is very small compared to its curvature 
radius then Eq. ^ implies that 
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V 

K 



0, 



(5) 



where v is the microscopic interface velocity and k is 
proportional to the mean curvature of the interface at 
each point. 



Assuming that the damping is homogeneous so that 
id = £d{t)i the macroscopic equation of motion may be 
written as 



where 



fv^pdV 



JpdV 



= 0. 



J v{l-v'^)KpdV 

(1 - v^)v j pdV 



(6) 



(7) 



p is the interface energy density and V is the physical 
volume. We have assumed that {v"^) = (the brackets 
denote a spatial average) in the derivation of Eq. (|6|. 
This assumption is valid in the relativistic regime up to 
first order in (1 — v) and it has a negligible impact in the 
nonrelativistic limit. Although it is possible to construct 
network configurations with the same v but quite differ- 
ent K, in most physically realistic situations it is sufficient 
to consider that R = k{v) [21] . 

Energy-momentum conservation in an expanding ho- 
mogeneous and isotropic background leads to 



^ + mip + p) = o, 



(8) 



where p = V^^ J pdV and p ~ V^^ J pdV are the av- 
erage pressure and density whose ratio is given by the 
equation of state parameter, w — p/p = —2/3 -I- v'^. Eq. 
([S]) can be generalized in order to account for additional 
energy loss mechanisms, due to interface collapse and 
friction: 



p + 3H{p + p) = - 



cv 
1 



If 



(9) 



Here, c is a phenomenological parameter which depends 
on the specific properties of the network and can be cali- 
brated using numerical simulations I19j , £f is the friction 
length scale, and L is the characteristic length of the 
interface network defined by p = cr/L, where a is the in- 
terface energy per unit area. Using Eq. ^ one obtains 



L 



L , 
Id 



cv . 



(10) 



Eqs. ^ and ( 10 1 constitute the VOS model for domain 
walls which was proposed in |22j and is here derived in 
detail. In particular the relation between the macroscopic 
parameter R and its microscopic counterpart k has been 
determined only in the present paper. 

Consider the case with = 00 and i? cx P (with 
a > 0). The attractor solution to the evolution Eqs. 



H 



6|) and ( 10 1 corresponds to a linear scaling solution with 
^ — t/a and 



L 

J 



I k{k + c) 
'ia{\ — a) 



(11) 



where k = kL is a dimensionless curvature parameter. 
A linear scaling solution can be attained only if a < 1 
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(-D > 1 or L < are unphysical). On the other hand, 
if a > 1/4 a Unear scahng solution is possible even with 
c = 0. If = oo then the interfaces are expected to be 
relativistic, as long as the expansion is sufficiently slow, 
except for very particular and unrealistic network con- 
figurations, with A; ^ 0. These results are supported by 
numerical simulations of cosmological domain wall net- 
works in a frictionless regime [19, 22 . On the other hand, 
if domain walls were to provide a significant contribution 
to the dark energy then including friction would not sig- 
nificantly help frustration due to the limited amount of 
energy with which domain walls can interact conserving 
energy and momentum |18l I19| . Hence, it is enough to 
consider the relativistic case with = oo in order to 
probe domain walls as a relevant dark energy candidate. 

In the context of phase field models the background is 
usually fixed (a = and H — 0) and the dynamics of the 
nonrelativistic interfaces is friction dominated. In this 
regime |w| <C \v/£f\ and Eq. ([s]) implies the microscopic 
relation v — nif. If the friction is homogeneous and time 
independent (ff = constant) then the interfaces move 
with a local velocity proportional to the local curvature 
(note that in the relativistic regime there is no direct 
relation between the local curvature and velocity). In 
the nonrelativistic case it is straightforward to find the 
following attractor solution to the evolution Eqs. ^ and 
pi: 



L = y/2?Mit , w=-i<l, (12) 

which is the usual t^^^ scaling for nonrelativistic inter- 
faces, implied by von Neumann' s law 

III. NUMERICAL RESULTS AND DISCUSSION 

A specific realization of a particular class of multicom- 
ponent models which, in the limit of an infinite number 
of coupled scalar fields, approaches the so-called ideal 
model, was developed and studied in detail [HI [19] in a 
cosmological context. It is characterized by a scalar field 
multiplet with iV^ real scalar fields and a scalar field po- 
tential with + 1 equivalent minima, equidistant in 
configuration space, 

V^H ^n^'-^oT with = 5^(0. (13) 

J=l i=l 

where the pi^ are the iV^ -I- 1 coordinates of the vacua of 
the potential, chosen to be the vertices of an {N + 1)- 
dimensional regular polyhedron with a distance between 
vertices equal to rg. The special properties of this class of 
models result from the fact that all domain walls have the 
same tension, so that only Y-type junctions are allowed. 
Also, the probability of fusion of adjacent domains as a 
result of domain wall collapse decreases proportionally to 
l/N^ for large N^, thus vanishing in the — oo limit. 



Although other classes of models allowing for X-type 
junctions |23] have also been investigated in Refs. [T91I24] 
it has been argued that the ideal model was the best 
model from the point of view of its potential to lead to 
frustration, a necessary condition for domain walls to 
provide a significant contribution to dark energy. On 
the other hand, in a non-relativistic friction-dominated 
regime and in the absence of expansion, the ideal model 
can be used to describe the dynamics of ideal grain 
growth. 

Here we compare the results of a set of high resolu- 
tion interface network simulations of the ideal class of 
models with = 10 (11 minima represented by 11 dif- 
ferent colors in Fig. 1), considering the relativistic and 
non-relativistic regimes studied in the previous section 
(we have taken a = 2/3 and a = 0, respectively). Fig. 1 
shows two 256^ snapshots taken from relativistic (left) 
and nonrelativistic (right) 1024^ two-dimensional inter- 
face simulations {d<f>/dz = 0) with periodic boundaries 
and identical initial conditions: At each grid point <f> was 
randomly assigned to one of the minima of the potential 
and (j) was initially set to zero (see Ref. [19j for a more de- 
tailed account of simulation details). The characteristic 
length of the network is roughly the same in both snap- 
shots (note that the corresponding physical times can be 
very different). Fig. 1 shows that although relativistic 
and nonrelativistic interface networks have very different 
dynamics, a single simulation snapshot does not clearly 
distinguish between the two regimes (we verified that this 
result is valid irrespective of the specific choice of A^^). 

This is shown more quantitatively in Fig. 2, which 
summarizes the results of ten relativistic and nonrela- 
tivistic snapshots with 1024 domains each, taken from 
1024^ numerical simulations. Again, the same initial 
conditions were considered for each pair of relativistic- 
nonrelativistic simulations. The blue stripe represents 
the interval /„ ± s„ where /„ is the fraction of the to- 
tal number of domains with n edges for either relativis- 
tic (upper panel) or nonrelativistic (lower panel) regimes 
and is the sample variance calculated using ten 1024^ 
domain wall network simulations. The red dots repre- 
sent the value of /„ in the nonrelativistic (upper panel) 
or relativistic cases (lower panel) for comparison. We 
verified that the average number of edges of an individ- 
ual domain, numerically calculated from the simulations, 
is equal to 6 in both cases, a characteristic of interface 
networks with Y-type junctions. Also, we observed no 
significant changes to /„ on subsequent snapshots. 

The distribution function /„ has been extensively used 
to characterize nonrelativistic two-dimensional patterns 
(see Ref. [1 for a review). Analytical predictions for 
/„ have been made in the context of soap froth and 
other similar systems [HI [HI [21]. In particular, the scale- 
invariant distribution function /„ for two-dimensional 
soap froths was found to be in very good agree- 
ment with the Monte Carlo simulations of the ideal Potts 
model |27| , corresponding to an infinite number of differ- 
ent domains. The most significant difference between the 
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FIG. 1: (Color online) Two 256^ simulation snapshots, with 
similar characteristic lengths, taken from 1024^ relativistic 
(left) and nonrelativistic (right) two-dimensional interface 
simulations of the ideal class of models with A'^^ — 10, with 
identical initial conditions. Different colors represent different 
minima. 



nonrelativistic simulation results of Fig. 2 (with = 10) 
and soap froth experimental results is that no two-sided 
domains are actuaUy observed in real soap froths (cor- 
responding to iV^ = oo). The fact that /„ is maximal 
for n = 5 is rather standard in relaxed two-dimensional 
systems with Y-type junctions, such as real soap froths. 

It was demonstrated in Ref. [KT that in order to pro- 
vide a significant contribution to the dark energy without 
generating exceedingly large CMB temperature fluctua- 
tions, a domain wall network must have a characteristic 
length L < lOkpc and a characteristic velocity v < 10~^ 
at the present time. That would be possible only if the 
relativistic distribution function in Fig. 2 was sharply 
peaked around n — 6. As shown in Fig. 2 this is clearly 
not the case. Hence, the presence of two-, three-, four-, 
and five-edge unstable domains is enough to prevent the 
frustration of the interface network. 

Fig. 2 shows that for n > 6 the number density of 
domains with a given number of edges is similar for rel- 
ativistic and nonrelativistic networks, in agreement with 
Fig. 1, although there are some differences for n < 6. In 
the relativistic case the velocity can never be much larger 
than V, since it is limited by the speed of light. On the 
other hand, in the nonrelativistic regime, v cx k£{ ^ v 
for highly curved interfaces, leading to the suppression 
of very small domains in the nonrelativistic regime. This 
can be seen in Fig. 1, and it is the main source of the 
differences between the relativistic and nonrelativistic re- 
sults shown in Fig. 2. Still, the similarity between rela- 
tivistic and nonrelativistic snapshots of simulations with 
identical initial conditions is striking, thus demonstrat- 
ing that it is the initial conditions rather than the specific 
dynamics of the cellular pattern which essentially deter- 
mine the large scale network properties at late times. 
Although in this paper only the results of interface simu- 
lations in two spatial dimensions are presented, we have 
verified that similar results also hold in the case of two- 
dimensional snapshots of three-dimensional interface net- 
work simulations. 
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FIG. 2: (Golor online) The blue (light gray) stripe represents 
the interval ± s„, where /„ is the fraction of the number 
of domains with n edges for either relativistic (upper panel) 
or nonrelativistic (lower panel) regimes and is the sample 
variance calculated using ten snapshots with 1024 domains 
each, taken from 1024^ numerical simulations of the ideal class 
of models with A^^ = 10. The red dots represent the value 
of /„ in the nonrelativistic (upper panel) or relativistic cases 
(lower panel) for comparison. 



These results suggest that a modified statistical version 
of von Neumann' s law might be valid for relativistic in- 
terface networks. This can be confirmed by determining 
the evolution of the comoving area of a single domain of 
n sides with physical time, d„ — ~R^^ § vdl, where dl is 
the comoving infinitesimal interface arc length (note that 
the physical length and area may be obtained by multi- 
plying the comoving ones by R and i?^, respectively). 
Ignoring the v term in Eq. ([s]) one has 



_ed_ 

i?2 



Krdl 




(14) 



where fee — Rk and 6i represent each of the n discon- 
tinuous angle changes at the vertices. If the interface 
network has only Y-type junctions which meet at an an- 
gle of 27r/3, then 



i?2 



27r 



27r 



3 i?2 



(n-6). (15) 



For scaling relativistic interface networks it is v and not v 
which is equal to zero and consequently von Neumann's 
law is valid only on average. The evolution of the total 
number of domains, Nd, with time can also be obtained 
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fromEq. (15). If /„ = then 



Nd 



(16) 



Here A — J^'^^i is the comovmg area of the entire 
system, A„ is the total comoving area occupied by do- 
mains with n edges, a — A/ No is the average comov- 
ing domain area, and a„ = (7„a is the average comoving 
area of a domain with n edges. If both /„ and are 
time independent, then Nd oc —i^N'jj/R^. Eqs. (15) 

J Y . —1/2 . . — ' 

and ( 16 ) imply that L cx RNj^ is proportional to t in 

the frictionless relativistic case cc H — a/t, with 
< a < 1) and proportional to t^/^ in the nonrela- 
tivistic friction dominated regime (^d — £i and R — 1), 



in agreement with Eqs. (11) and (12 1. Therefore, the 



main difference between relativistic and nonrelativistic 
thin featureless interface networks with the same C and 
initial conditions is associated with the scalings of the 
macroscopic parameters v and L with physical time, i, 
rather than with the properties of the cellular patterns 
taken from single simulation snapshots with a fixed L. 
Given that frustration, a necessary condition for a do- 
main wall network to be able to accelerate the expansion 
of the universe, cannot happen if there are unstable do- 
mains with less than six edges, laboratory tests, such as 
those probing the dynamics of real foam coarsening [7] 
(see also Ref. [2] and references therein) can be used to 



rule out domain wall networks as a relevant dark energy 
source. 



IV. CONCLUSION 

In this paper we have provided a rigorous derivation 
of the VOS domain wall model, showing that it may be 
used to successfully describe the dynamics of both rela- 
tivistic and nonrelativistic interfaces. Although the evo- 
lution of the characteristic length and velocity of inter- 
face networks with physical time is very different in these 
regimes, we have shown that a single simulation snap- 
shot does not clearly distinguish relativistic and nonrela- 
tivistic networks with the same characteristic length, the 
reason being that a modified version of von Neumann's 
law also applies (in a statistical sense) to the relativistic 
regime. We have argued that this allows the use of lab- 
oratory tests with nonrelativistic interfaces to rule out 
cosmological domain wall networks as a significant dark 
energy source. 
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